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Populations of hematopoietic stem cells and progenitors are quite heterogeneous and consist of multiple cell 
subsets with distinct phenotypic and functional characteristics. Some of these subsets also appear to be 
interconvertible and oscillate between functionally distinct states. The multipotent hematopoietic cell line 
EML has emerged as a unique model to study the heterogeneity and interconvertibility of multipotent 
hematopoietic cells. Here we describe extensive phenotypic and functional heterogeneity of EML cells which 
stems from the coexistence of multiple cell subsets. Each of these subsets is phenotypically and functionally 
heterogeneous, and displays distinct multilineage differentiation potential, cell cycle profile, proliferation 
kinetics, and expression pattern of HSC markers and some of the key lineage-associated transcription 
factors. Analysis of their maintenance revealed that on a population level all EML cell subsets exhibit 
cell-autonomous interconvertible properties, with the capacity to generate all other subsets and re-establish 
complete parental EML cell population. Moreover, all EML cell subsets generated during multiple cell 
generations maintain their distinct phenotypic and functional signatures and interconvertible properties. 
The model of EML cell line suggests that interconvertible multipotent hematopoietic cell subsets coexist in a 
homeostatically maintained dynamic equilibrium which is regulated by currently unknown cell-intrinsic 
mechanisms. 



Numerous studies have established phenotypic and functional heterogeneity within populations of embry- 
onic stem (ES) cells, adult neural, intestinal and hematopoietic stem cells (HSCs) and cancer stem cells 
(CSCs), which arises from the coexistence of different stem cell subsets 118 . Moreover, a number of studies 
reported the existence of interconvertible cell subsets among ES cells, intestinal stem cells and HSCs, which 
oscillate between several metastable states, as well as bidirectional interconversions between mammary and 
intestinal CSCs and non-CSCs 19 " 28 . These findings reformed our perception of stem cells as a functionally 
uniform pool to that of a dynamic pool of multiple stem cell subsets. They are also resurrecting the question 
are stem cells distinct and permanent cellular entities or do they exist in several functional sates with distinct 
phenotypic and functional features 28,29 ? 

The increasingly better characterized heterogeneity among HSCs supports the notion that the HSC pool 
consists of different cell subsets with distinct phenotypic and functional characteristics 517 . In addition, new 
studies reported further phenotypic and functional heterogeneity among early and lineage-restricted hemato- 
poietic progenitors 14,15 . The range of expression of certain markers (e.g. CD34, CD150, CD229, CD41) on HSCs 
and differential Hoechst 33342 efflux capacities enabled the identification of functionally distinct HSC subsets, 
which differ in their self-renewal potential, life-span, cycling status and differentiation potential. More impor- 
tantly, it appears that some of these HSC subsets are interconvertible and oscillate between functionally distinct 
states 8,10,11,19,22,23 . 

Studies by the Ogawa group were the first to reveal that mouse long-term repopulating HSCs (LTR-HSCs) 
oscillate in vivo between CD34~ and CD34 +/low states 18,19 . Other studies suggested that mouse LTR-HSCs oscillate 
between dormant CD34~ and activated CD34 + states 22 . Based on differential expression of CD150 marker and 
Hoechst dye efflux capacity, several studies have shown that HSCs can be divided into two functionally distinct 
and potentially interconvertible lineage-biased HSC subsets: (1) a CD150 hl subset with high propensity to 
generate myeloid progeny, and (2) a CD150 low subset with robust potential for differentiation into lymphoid 
lineages 6 - 8,13,14,23 . 
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Consequently, HSCs are now perceived as a complex and dynamic 
pool of functionally heterogeneous and sometimes interconvertible 
cell subsets, rather than being a functionally uniform cell popu- 
lation 6 " -1016,23 . However, the characteristics and regulation of inter- 
convertible HSC subsets remain to be fully explored. Studying HSC 
interconvertibility in vivo or in vitro is complicated due to difficulties 
in capturing oscillating cell subsets in vivo and long-term in vitro 
maintenance of functional HSCs 4 . 

The murine multipotent hematopoietic cell line EML has emerged 
as a unique model to study heterogeneity and interconvertibility of 
multipotent hematopoietic cells 30 " 34 . The EML cell line is SCF- 
dependent and was derived by retroviral expression of truncated, 
dominant-negative form of human retinoic acid receptor 
RARa403 in BM cells from 5-fluorouracil treated BDF1 mice 30 . 
Importantly, EML cell line has a well-documented multilineage 
(erythroid, myeloid, lymphoid) differentiation capacity, and is an 
established in vitro surrogate for multipotent hematopoietic 
cells 30,35 " 41 . In the presence of SCF EML cells undergo proliferative 
self-renewal, while in response to specific cytokines and in presence 
of specific stromal cells they differentiate into cells of the erythroid, 
granulocyte-macrophage, megakaryocytic, B cell and T cell lineages 
(Supplementary Fig. SI) 30,35 " 41 . 

Multiple studies have drawn parallels between EML cells and adult 
HSCs and multipotent progenitors (MPPs) due to their phenotypic, 
functional and molecular features. For that reason EML cell line is 
being increasingly used as an important and biologically relevant 
model to study the regulation of multipotent hematopoietic cell 
self-renewal, lineage commitment and differentiation 30,35 " 41 . 

Previous studies have shown that EML cell line is phenotypically 
and functionally heterogeneous 30 " 32,34,38,39,41 . Ye et al. reported that 
EML cell line consists of Lin"CD34 + replicating cells and 
Lin"CD34" nonreplicating cells which do not grow in the presence 
of SCF 31 . Furthermore, that study reported that the CD34 + cell popu- 
lation generates a CD34" nonreplicating population, and maintains 
both cell populations within EML cell line 31 . The subsequent study by 
the same group expanded their observations that Lin"CD34 + cell 
subset grows rapidly in SCF and can regenerate entire population 
of EML cells, whereas Lin"CD34" subset grows slowly if at all in SCF 
but can be propagated in the presence of IL-3 33 . 

Chang et al. have shown that EML cells exhibit ~ 1,000-fold range 
in the expression level of HSC marker Sca-1 32 . By purification and 
culture of EML cells with the highest, middle and lowest —15% Sea- 1 
expression level (referred to as Sca-1 1 "™, Sca-l mld and Sca-l hl8h sub- 
sets), that study has shown that on a population level each Sca-1 
subset is interconvertible and regenerates the parental EML cell 
population and its phenotypic and functional traits, including the 
expression pattern of Sca-1 marker, differentiation rates and gene 
expression profiles. In addition, using methylcellulose cloning Chang 
et al. generated random single cell-derived clones from total EML 
cells that expressed varying levels of Sca-1, which over a period of 8 
weeks regenerated the original broad expression profile of Sca-1, 
typical for EML cell line. Based on these findings Chang et al. pro- 
posed that interconvertible states of EML cells reflect slowly fluctuat- 
ing transcriptomes that may govern reversible priming and cell fate 
decision 32 . 

A study by Pina et al. characterized interconvertible properties and 
lineage commitment status of Sca-1 1 " EML cells (encompassing the 
15% lowest-staining cells) and Sca-l hl EML cells (encompassing the 
15% highest-staining cells) 34 . Although on a population level both 
Sca-1 1 " and Sca-l hl EML cells have culture-reconstituting capacity, 
subsequent analysis of the CD34 + and CD34" fractions of Sca-l'° 
EML cells revealed that Sea- 1'° population consists of erythroid-com- 
mitted Seal 1 " CD34" cells with no culture-reconstitution capacity 
and a culture-reconstituting Seal 1 " CD34 + cells with both erythroid 
and myeloid differentiation potential. By analyzing the so called 
culture-reconstituting potential of single Sca-1 1 ", Seal 1 " CD34" and 



Seal 1 " CD34 + cells over a period of several days, Pina et al. reported 
that 17% of Sca-1 1 ", 2.5% of Seal 1 " CD34" and 35% of Seal 1 " CD34 + 
cells generated more than 100 cells, which reconstituted broad Sca-1 
expression. These findings led to the conclusion that Sca-1 subsets of 
EML cells reflect oscillation between self-renewing and committed 
cell states 34 . 

To gain further insight into the heterogeneity and interconvert- 
ibility of multipotent hematopoietic cells, in this study we used the 
broad expression range of the HSC marker CD34 on EML cells to 
analyze in more depth (1) the phenotypic and functional heterogen- 
eity within EML cell line, and (2) the characteristics and maintenance 
of interconvertible EML cell subsets. 

Results 

EML cell line is phenotypically very heterogeneous and consists of 
multiple cell subsets with distinct expression pattern of HSC 
markers. To gain further insight into heterogeneity within EML 
cell line we examined the expression of a panel of HSC cell surface 
markers 5,42 " 49 . Flow cytometry analysis has shown that all or nearly all 
EML cells are Sca-1 + c-kit + Flk2" and alpha2-integrin (CD49b)", but 
exhibit heterogeneous expression of CD34, CD27, VEGFR1, CD201 
and CD150 (Fig. la,b and Supplementary Fig. S2). 

Division of HSCs into CD150 neg , CD150 1 "", CD150 inl and 
CD150 hlgh subsets based on the range of CD150 marker expression 
has revealed that these subsets display distinct in vivo functional 
features and differential lineage transcriptional priming 6,15 . Further- 
more, expression range of CD34 marker is used to distinguish LTR- 
HSCs from short-term repopulating HSCs (STR-HSCs), and perhaps 
dormant from activated LTR-HSCs 22,42,44 . Similarly, fractionation of 
EML cells into Sca-l'° w , Sca-l m,d and Sca-l h,gh subsets based on the 
broad expression range of Sea- 1 marker has enabled characterization 
of EML Sca-1 cell subsets 31,32,34 . 

EML cells also exhibit a broad expression range of CD34, similar to 
that of Sca-1 (Fig. la,c) 32 . Using equal increments of fluorescence 
intensity we divided EML cells into CD34 neg , CD34'° W , CD34 int and 
CD34 high cell subsets (referred to in further text as EML CD34 cell 
subsets), designating the cells that do not express CD34, and 
cells expressing low, intermediate and high levels of CD34 
(Fig. lc). This approach enabled us to reproducibly and consistently 
purify these EML cell subsets, and study their heterogeneity and 
interconvertibility. 

Although the frequency of each EML CD34 cell subset fluctuates 
slightly, their average frequency remains quite stable throughout 
multiple cell passages. Further analysis of EML CD34 cell subsets 
revealed phenotypic heterogeneity and distinct expression profiles of 
CD27, CD201, VEGFR1 and CD150 markers within each cell subset. 
For example, CD34 neg and CD34 low subsets contain the highest per- 
centage of VEGFR1 + cells and the lowest percentage of CD150 + 
and CD27 + cells, whereas CD34 int and CD34 high subsets contain 
the highest percentage of CD150 + and CD201 + cells (Fig. Id and 
Supplementary Fig. S3). 

Flow cytometry analysis of co-expression of CD34 and Sca-1 mar- 
kers revealed further differences between EML CD34 subsets. The 
CD34 neg subset reproducibly exhibits the broadest range of Sca-1 
expression, which gradually narrows from CD34 low to CD34 mt and 
CD34 hlgh cell subsets (Supplementary Fig. S4). Although the majority 
of cells in the CD34"' e subset are Sca-1 + , this is the only subset that 
also contains Sca-l neg cells. In contrast, the cells in CD34 hlgh subset 
express the highest levels of Sca-1 (Supplementary Fig. S4). 

EML CD34 cell subsets are functionally diverse and display distinct 
multilineage differentiation patterns, cell cycle profiles and 
proliferation kinetics. To investigate the functional characteristics 
of EML CD34 cell subsets, we analyzed their multilineage differen- 
tiation potential, cell cycle profile and proliferation kinetics. 
Multilineage differentiation of FACS purified CD34 neg , CD34 low , 
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Figure 1 | The multipotent hematopoietic cell line EML is phenotypically heterogeneous and consists of several cell subsets with distinct expression 
patterns of HSC markers, (a) Representative flow cytometry analysis of the expression of a panel of HSC markers on EML cells showing that almost 
all EML cells are Sca-1 + c-kit + Flk2~ and alpha2-integrin (CD49b)", but exhibit heterogeneous expression of CD34, CD27, VEGFR1, CD201 and CD150 
(« = 6). (b) Average percent of EML cells expressing CD27, CD201, VEGFR1, CD150and CD34 markers (n = 6). Data are mean ± SD. (c) Representative 
flow cytometry analysis showing division of EML cells into CD34 neg , CD34 low , CD34 tal and CD34 high cell subsets based on the broad range of 
expression of the CD34 marker (n = 15). (d) Histograms show the results of flow cytometry analysis of CD27, CD201, VEGFR1 and CD 150 expression on 
CD34 nes , CD34 low , CD34 int and CD34 hl8h EML cell subsets (n = 3-9). Data are mean ± SD. * P value < 0.05, ** P value < 0.001. 



CD34' nt and CD34 hlfih cell subsets (-99% purity, >95% viability) was 
assessed by colony- forming cell (CFC) assays for hematopoietic 
progenitors 37,38 . CFC assays revealed that only a very small percen- 
tage of CD34"' e cells forms erythroid (BFU-E) and granulocyte- 



macrophage (CFU-GM) colonies (Fig. 2a). The CD34 nee cells that 
do not form colonies are visible as viable single cells for up to 5 days 
after plating. During that period these cells do not divide at all, and 
around day 5 they begin to die and are visible as cell debris. This 
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Figure 2 | EML CD34 cell subsets exhibit distinct multilineage differentiation patterns, cell cycle profiles and proliferation kinetics, (a) Colony- 
forming cell (CFC) assay analysis of the differentiation capacity of total EML cells and purified CD34 ne6 , CD34 low , CD34 int and CD34 hl8h EML cell subsets 
into erythroid, megakaryocytic and granulocyte-macrophage lineages (n = 15). The histograms and the table show the numbers of granulocyte- 
macrophage (CFU-GM), megakaryocytic (CFU-Meg) and erythroid (BFU-E) colonies generated by total EML cells and each EML CD34 cell subset 
(10 3 cells/well) in the presence of cytokines IL-3, Tpo andEpo, respectively. Data are mean ± SD. * P value < 0.05, ** P value < 0.001. (b) Analysis of cell 
cycle profiles of total EML cells and CD34 ne8 , CD34 low , CD34 int and CD34 hish EML cell subsets using Propidium Iodide, Hoechst 33342 and Ki-67 vs 7-AAD 
staining (n = 3-5). Data are shown as mean ± SD. * P value < 0.05, ** P value < 0.001. (c) Flow cytometry analysis of proliferation kinetics of total EML 
cells by CFSE labeling. The grey-shaded histogram on the left side of the panel shows unlabeled cell auto-fluorescence and the limits of detectable cell 
divisions. The histogram outlined in red on the right side of the panel shows CFSE signal in EML cells immediately after the labeling. The overlay of CFSE 
signal at 24 hour time points following labeling shows EML cell division history. During each cell division CFSE is partitioned equally among daughter 
cells, resulting in a corresponding uniform decrease in CFSE fluorescence among all dividing cells. Synchronous cell divisions result in a corresponding 
uniform decrease in CFSE signal among all dividing cells and narrow fluorescence histograms. (n = 2). (d) Analysis of proliferation kinetics of CD34 neg , 
CD34 low , CD34'" 1 and CD34 h ' 8h EML cell subsets by CFSE labeling. The overlay of CFSE signal at 24 hour time points following labeling shows division 
history and proliferation kinetics of each EML cell subset. The cells within CD34" es subset exhibited more asynchronous and slower division at each time 
point (see also Supplementary Fig. S5), as shown by varying levels or CFSE fluorescence resulting in broader and overlapping histograms at different time 
points after CFSE labeling (n = 4). 
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observation indicates that cytokines IL-3, Epo or Tpo do not induce 
differentiation or support survival of the majority of cells within the 
CD34 neg subset. In contrast, the percentage of cells forming 
granulocyte-macrophage, erythroid and megakaryocyte (CFU- 
Meg) colonies increases from CD34 low to CD34 hlBh cells (Fig. 2a), 
indicating an increased content of colony-forming progenitors in 
these EML cell subsets. Notably, these findings are very similar to 
a report that mouse Lin" Sca-1 + c-kit + (LSK) CD34 high BM cells 
differentiate efficientiy into myelo-erythroid lineages, whereas 
CD34 low and CD34~ BM cells differentiate minimally or not at all 19 . 

Analyses of the cell cycle profiles by PI, Hoechst 33342, and Ki-67/ 
7-AAD staining (Fig. 2b) have all shown that each EML CD34 cell 
subset contains a significantly different percentage of cells in G0/G1 
and S/G2/M phases when compared to other subsets. While a large 
proportion of CD34 neg cells are in G0/G1 phase, the percentage of 
cells in S and G2/M gradually, but significantly, increases from 
CD34 iow tQ CD3 4i.igi. ce ii SUDset ( Fig . 2 b). 

To examine their proliferation kinetics, total EML cells and puri- 
fied EML CD34 cell subsets were labeled with CFSE and their divi- 
sion evaluated at 24 hour time points for 5 days (Fig. 2c,d). The CFSE 
assay revealed that total EML cells divide quite synchronously, with 
only a small proportion of cells dividing more slowly (Fig. 2c and 
Supplementary Fig. S5). Remarkably, CFSE labeling uncovered the 
existence of two distinct cell populations within total EML cells: a 
prominent large CFSE Mam population with uniformly high level of 
CFSE signal, and a much smaller CFSE Sate "" e population with lower 
level of CFSE signal immediately after labeling (Supplementary Fig. 
S5b). The CFSE Mam population is very similar to total EML cells in 
terms of their CD34 expression pattern and the relative frequency of 
CD34 ne8 , CD34 low , CD34 tol and CD34 high cell subsets. In contrast, the 
CFSE sa.emte p opu i ation 0 f EML ce u s consists mainly of CD34 neg and 
CD34 low cells (Supplementary Fig. S5b,c). 

The cells within the CD34 hlBh and CD34 int subsets also divide quite 
synchronously, whereas a small fraction of cells in CD34 low subset 
divides asynchronously (Fig. 2d and Supplementary Fig. S6a). In 
contrast, the cells within the CD34 neg subset divide more asynchro- 
nously, with a portion of cells dividing slower than the rest (Fig. 2d 
and Supplementary Fig. S6a). 

CFSE labeling of EML CD34 subsets also uncovered the existence 
of distinct CFSE Main and CPSE Satellit ' ! subpopulations mainly within 
the CD34 neg and CD34 low cell subsets (Supplementary Fig. S6b). The 
frequency of CFSE Satemte cells is the highest in the CD34 n,!g subset 
(22.68 ± 6.7%), while the CD34 low cell subset contains lower fre- 
quency (9.4 ± 5.01%) of CFSE Salemte cells. On the other hand, the 
CD34 mt and CD34 hlgh cell subsets contain a very small percentage of 
CFSE sa.diite ceUs q 03 ± 0.56% and 0.45 ± 0.37%, respectively) 
(Supplementary Fig. S6b). Temporal tracking of the CFSE signal in 
the progeny of each EML CD34 cell subset in relation to CFSE Mam and 
CFSE Sa,<:lllte populations further illustrates asynchronous and slower 
proliferation of cells in the CD34™ g and CD34 low subsets in compar- 
ison to CD34 tot and CD34 hlgh cell subsets (Supplementary Fig. S6c). 

Overall, CFSE assays revealed additional heterogeneity within 
EML cell line, and additional differences between the CD34 neg / 
CD34 low and CD34 in 7CD34 high cell subsets, which may be due to 
(1) increased CFSE efflux capacity of CFSE SaMU " e cells, or (2) reduced 
activity or amount of intracellular esterases required to remove acet- 
ate groups and render CFSE fluorescent. 

Cumulatively, these studies revealed additional phenotypic and 
functional diversity within EML cell line, and support the notion 
that in terms of their phenotypic and functional features the 
CD34 neg and CD34 low subsets are most similar to each other and share 
many phenotypic and functional features, whereas on the other end 
the CD34 tal and CD34 high subsets share most similarities. 

Hematopoietic transcriptional signatures of EML CD34 cell sub- 
sets. To start identifying molecular differences between EML CD34 



cell subsets we analyzed transcriptomes in (1) three independent 
samples of total EML cells, and (2) three independently purified 
samples of CD34 neg cells and combined CD34 int and CD34 Mgh 
subsets, referred to in text as CD34 + cells. To compare hematopoie- 
tic transcriptional signatures of these EML cell subsets we specifically 
examined expression of (1) genes important for overall hemato- 
poiesis and HSC maintenance and function, (2) hematopoietic cell 
surface markers, receptors and adhesion molecules, and (3) genes 
associated with erythroid, megakaryocyte, myeloid and lymphoid 
lineage development and specification 50 52 (Supplementary Tables 
1,2,3). 

Analysis of microarray data detected little or no difference in 
expression (£ 1.5-fold change) of the majority of genes examined 
in CD34 neg and CD34 + EML cells (Supplementary Tables 1,2,3). 
Genes whose expression was higher in CD34 neg EML cells included: 
GATA-1, junctional adhesion molecule 4 (JAM4), CD38 (another 
mouse HSC marker), Tie2/TEK, HSC-derived tyrosine kinase STK 
(RON), a receptor for Macrophage-stimulating protein (MSP), Evil, 
Hesl, Klfl (Kruppel-like factor 1 erythroid) and VWF 
(Supplementary Tables 1,2,3). Genes found to be down-regulated 
in CD34 neg EML cells included: PU.l, C/EBPoc, CD34, CD27 
(Fig. lc, d), JAM3, Idl, Tnfrsfloc, Lck, TdT and Tcf7 (Supple- 
mentary Tables 1,2,3). The expression of transcript for Tcf7 gene 
(a transcriptional activator that plays an important role in lympho- 
cyte differentiation) was recently found to be significantly down- 
regulated in Sca-1~ CD34~ EML cells in comparison to Sca-1 + 
CD34 + cells 38 . 

Of particular interest among the genes differentially expressed 
between CD34 neg and CD34 + EML cells are transcription factors 
GATA-1, PU.l and C/EBPa 53 57 (Fig. 3). PU.l and GATA-1 play 
key roles in myeloid and erythroid lineage specification and restric- 
tion 53,54 . Their reciprocal repression controls the erythroid versus 
myeloid fate decisions 53,54 . C/EBPoc controls developmental fate of 
multipotent progenitors and HSCs by actively promoting myeloid 
differentiation and inhibiting erythroid differentiation 55,56 . The latest 
study reported that C/EBPoc is required for myeloid lineage priming 
of HSCs, and that deletion of C/EBPoc results in loss of self- renewal 
and HSC exhaustion 57 . 

Previous studies reported that Sca-l low and Sca-l low CD34" EML 
cells express significantly higher levels of GATA-1 transcript and 
protein than Sca-l med , Sca-l high , and Scal low CD34 + EML cells 32,34 . 
Chang et al. have also observed that Sea- l h,gh and Sea- l med EML cells 
express significantly higher levels of PU.l transcript and protein than 
Sca-l low EML cells 32 . Earlier microarray studies also revealed much 
higher expression of C/EBPoc transcript in Seal + CD34 + EML cells in 
contrast to Scar CD34" EML cells 41 . 

Our microarray data indicated that the expression of GATA- 1 is 
higher in CD34 neg EML cells and the expression of PU.l and C/EBPoc 
is higher in CD34 + cells (Supplementary Tables 1,2,3). Thus, we first 
analyzed the expression of GATA-1, PU.l, C/EBPoc genes in purified 
CD34 neg , CD34 low , CD34 lnt and CD34 high EML cell subsets by quant- 
itative real-time RT-PCR (qRT-PCR) (Fig. 3a). 

qRT-PCR analysis revealed significant differences in the express- 
ion of GATA- 1 and PU. 1 transcripts among EML CD34 subsets with 
(1) the expression of GATA-1 being significantly higher in the 
CD34 neg and CD34 low cells than in the CD34"" and CD34 high cells, 
and (2) the expression of PU.l being the highest in the CD34 h,gh cells 
and significantly lower in the CD34" 11 , CD34 10 ™ and CD34 neg subsets 
(Fig. 3a). qRT-PCR has also shown that CD34 Mgh and CD34 tot cells 
express similar but significantly higher levels of C/EBPoc transcript 
than CD34 neg and CD34 low cells (Fig. 3a). 

The microarray data also indicated that CD34 neg EML cells express 
higher level of junctional adhesion molecule 4 (JAM4), which is 
transcribed in LSK CD34~ and LSK CD34 + BM cells containing 
HSCs and MPPs, but not in more differentiated BM cells 58 
(Supplementary Tables 1,2,3). qRT-PCR has shown that the express- 
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Figure 3 | Analysis of hematopoietic transcriptional signatures in EML CD 34 cell subsets. To identify molecular differences between EML CD34 cell 
subsets the transcriptomes of purified CD34 ne6 cells and combined CD34 1 " 1 and CD34 h ' gh subsets, refrerred to in the text as CD34 + cells, were analyzed 
using Agilent microarrays. This analysis focused specifically on expression of ( 1 ) genes important for overall hematopoiesis and HSC maintenance 
and function, (2) hematopoietic cell surface markers, receptors and adhesion molecules, and (3) genes associated with erythroid, megakaryocytic, 
myeloid and lymphoid lineage specification and development, (a) The quantitative real-time RT-PCR expression analysis of transcription factors GATA- 
1, PU. 1 and C/EBPa, and adhesion molecule JAM4 in purified CD34™ 8 , CD34 low , CD34 int and CD34 hi8h EML cells {n = 3). Data are shown as mean + SD. * 
P value < 0.05, ** P value < 0.01. (b) Representative Western analysis of GATA-1, PU.l, C/EBPa and JAM4 protein expression in total EML cells and 
purified CD34 ne6 , CD34 low , CD34 int and CD34 hl8h EML cell subsets. Cell lysates from U-937 cells were used as positive controls (n = 3). The histograms 
show the levels of GATA- 1 , PU. 1 , C/EBPa and JAM4 proteins in each EML CD34 cell subset normalized to GAPDH protein levels. The levels of GATA- 1 , 
PU.l, C/EBPa, JAM4 and GAPDH proteins in each sample were measured by densitometric scanning and analyzed with ImageJ software. The 
staining of gels with Coomassie Brilliant Blue shows loading of proteins from total EML cells and EML CD34 cell subsets. 



ion of JAM4 transcript is indeed significantly higher in CD34 neg cells 
than in CD34 low , CD34 1 "' and CD34 hish cells, thus confirming micro- 
array data (Fig. 3a). 

The expression of GATA-1, PU.l, C/EBPa and JAM4 in purified 
EML CD34 cell subsets was further examined by Western analysis. 
The expression of GATA-1 protein was found to be higher in the 
CD34 ne s and CD34 low cells than in the CD34 tot and CD34 h ' 8h cells, 
while the expression of PU. 1 protein was found to be higher in the 
CD 3 4 high and CD34 tot ce u s trian m me CD34 low and CD34 ne e cells 

(Fig. 3b). 

Notably, although the CD34 h ' 6h and CD34 int cells express much 
higher levels of C/EBPa transcript than CD34 ne6 and CD34 low cells, 
Western analysis detected no significant difference in the expression 
of C/EBPa protein among EML CD34 cell subsets (Fig. 3a,b). 

On the other hand, the expression pattern of JAM4 protein in 
purified EML CD34 cell subsets was very similar to expression pat- 
tern of JAM4 transcript, with the expression of JAM4 protein being 
the highest in the CD34 nt ' s cells and the lowest in the CD34 h ' 8h cells 
(Fig. 3a,b). 

EML cell line consists of interconvertible cell subsets which are 
maintained and coexist in a dynamic but stable equilibrium. Using 
phenotypic and functional signatures of EML cell subsets as a blue- 
print, we interrogated the developmental relationships and main- 
tenance of all EML CD34 cell subsets during multiple passages and 
cell generations. Multiple batches of each EML CD34 cell subset were 
purified and subsequently cultured independently (Fig. 4a). To 
minimize the cross-contamination between cell subsets and vari- 
able cell viability, only the cells from sorts yielding the highest cell 



purity (>99%) and viability (>95%) were used in these experiments. 
Multiple independent cultures of purified EML CD34 cell subsets 
were maintained by serial cell passage for up to 14 days, and 
expression of CD34 on the progeny of each subset analyzed every 
2 days (Fig. 4a). 

The temporal analysis of the primary progeny of each EML cell 
subset has shown that already by day 12 of culture each EML CD34 
cell subset gave rise to all other cell subsets and re-established the 
broad expression pattern of CD34 marker typical for parental EML 
cell population (Fig. 4a and Supplementary Fig. S7). These results are 
quite similar to previously reported regeneration of the parental EML 
cell population and its broad range of Seal marker expression by 
purified Sca-1 cell subsets 31 . However, the temporal kinetics with 
which each EML CD34 cell subset generated other subsets and re- 
established the parental EML cell population differed, with CD34 neB 
cells exhibiting the slowest and CD34 h,gh cells showing the fastest 
kinetics (Fig. 4a and Supplementary Fig. S7). 

In another set of experiments we analyzed expression of HSC 
markers on and cell cycle profiles of total day 14 progeny of inde- 
pendently cultured parental EML cells and purified CD34 cell sub- 
sets. The expression patterns of CD27, VEGFR1, CD201 and CD150 
markers on total progeny of parental EML cells and each CD34 cell 
subset were very similar to and not significantly different from the 
pattern on parental EML cells (Fig. 4b and lb). Moreover, all new 
CD34 cell subsets, generated either by parental EML cells or by each 
parental CD34 cell subset, exhibited the heterogeneous expression 
patterns of CD27, VEGFR1, CD201 and CD150 markers, typical of 
the parental EML CD34 cell subsets (Fig. 5a,b and Fig. Id). For 
example, both among parental EML cell subsets and their progeny, 
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Figure 4 | Temporal analyses of the progeny of each EML CD34 cell subset showing that each EML cell subset has the capacity to generate all other 
CD34 cell subsets and re-establish the entire parental EML cell population, (a) Multiple batches of CD34 nes , CD34 low , CD34'" 1 and CD34 hl6h EML cell 
subsets were purified by FACS. To minimize the cross-contamination between cell subsets and variable cell viability, only the cells from sorts yielding the 
highest cell purity (>99%) and viability (>95%) were used in these experiments. Independent cultures of each purified CD34 cell subset were maintained 
by serial cell passage for 14 days. Representative flow cytometry density plots showing temporal changes of the CD34 expression on the progeny of each 
EML CD34 cell subset over 14 days of culture, and re-establishment of the broad expression range of CD34 marker typical for parental EML cell 
population ( n = 3). (b) The parental EML cells and purified CD34 ne8 , CD34 low , 0034"" and CD34 M8h cell subsets were cultured independently for 14 days, 
and the expression of CD27, VEGFR1, CD201 and CD 150 markers on their total day 1 4 progeny analyzed by flow cytometry (n = 2) . Data are mean ± SD. 
(c) The cell cycle profiles of day 14 total progeny of parental EML cells and each EML CD34 cell subset were analyzed by flow cytometry using Propidium 
Iodide {n = 2). Data are mean ± SD. 
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Figure 5 | Progeny of total EML cells and purified EML CD34 cell subsets consists of new EML CD34 cell subsets that recapitulate the parental 
expression pattern of HSC markers. Multiple independent cultures of total EML cells and purified parental CD34 nes , CD34 low , CD34' nt and CD34 h,gh cell 
subsets were maintained by serial cell passage for 14 days, and expression of several HSC markers on their day 14 progeny analyzed by flow cytometry. 
Total progeny of parental EML cells and each parental EML CD34 cell subset consisted of new CD34™ 8 , CD34 low , 0034"' and CD34 high cell subsets, (a) The 
histogram represents the results of flow cytometry analysis of CD27, VEGFR1, CD201 and CD150 HSC marker expression on new CD34 nes , CD34 low , 
CD34"" and CD34 1,18h cell subsets generated by parental EML cells after 14 days of culture (n = 2). Data are mean ± SD. (b) The histograms represent the 
results of flow cytometry analysis of CD27, VEGFR1, CD201 and CD150 HSC marker expression on new CD34 neB , CD34 low , CD34' nt and CD34 hi8h cell 
subsets generated by each of the parental EML cell susbets after 14 days of culture (n = 2). Data are mean ± SD. 



the CD34 neB and CD34 low cells in general contained higher percen- 
tages of VEGFR1 + cells and lower percentages of CD150 + cells than 
CD34 in ' and CD34 high cells, while the expression pattern for CD27 
and CD201 is opposite (Fig. 5a,b and Fig. Id). 

Although the overall pattern of HSC marker expression is pre- 
served, the mean percentage of cells expressing CD27 and CD201 is 
much more homogeneous among the progeny of each subset than 
among parental cell subsets. Also, the CD34 nee progeny of each par- 
ental subset contains higher percentage of VEGFR1 + cells than the 
parental CD34 neg cell subset (Fig. 5b and Fig. Id). These differences 



could stem from our observations (1) that the percentage of cells 
within each EML subset that express these markers fluctuates 
throughout culture (Fig. Id), and (2) that it could take up to 3 weeks 
for the cumulative progeny of each CD34 cell subset to reach stable 
expression of these markers which is by that time very similar to 
typical expression of HSC markers on parental cell subsets. 

Equally important, analysis of cell cycle profile of day 14 total 
progeny of the parental EML cells and each EML CD34 cell subset 
by PI staining has also shown (1) that the cell cycle profiles of the 
cumulative progeny of parental EML cells and each CD34 cell subset 
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Figure 6 | Progeny of total EML cells and purified EML CD34 cell subsets consists of new EML CD34 cell subsets that recapitulate the parental cell 
cycle profiles. Multiple independent cultures of total EML cells and purified parental CD34 ne8 , CD34 low , CD34 ml and CD34 hl8h cell subsets were maintained 
by serial cell passage for 14 days, and cell cycle profiles of their day 14 progeny were analyzed by flow cytometry using Propidium Iodide. Total 
progeny of parental EML cells and each parental EML cell subset consisted of new CD34 ne6 , CD34'° W , CD34 mt and CD34 hl8h cell subsets, (a) The histogram 
depicts cell cycle profile analysis of new CD34 ne6 , CD34 low , CD34 lnt and CD34 h ' 6h cell subsets generated by parental EML cells after 14 days of culture 
(n = 2). Data are mean ± SD. (b) The histograms depict cell cycle profile analysis ofnewCD34 neg , CD34 low , CD34"" and CD34 1,,gh cell subsets, generated by 
each parental EML CD34 cell subset after 14 days of culture (n = 2). Data are mean ± SD. 



are very similar and closely resemble cell cycle profiles of the parental 
population of EML cells (Fig. 4c and Fig. 2b), and (2) that the cell 
cycle profiles of new CD34 n<:g , CD34 low , CD34 int and CD34 hlBh cell 
subsets, generated either by parental EML cells or by each parental 
CD34 cell subset, are very similar and closely resemble the cell cycle 
profiles of the parental CD34 cell subsets (Fig. 6a,b and Fig. 2b). 



The interconvertible characteristics of each EML CD34 cell subset 
were further tested and analyzed through serial cell sorting, culture 
and analysis of their successive progeny over several generations. 
EML cells were sorted into parental CD34 nee , CD34 low , CD34"" and 
CD34 hieh subsets (cell purity > 99% and viability > 95%) which were 
cultured independently for 14 days (Fig. 7a). Analysis of CD34 
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Figure 7 | Serially generated EML cell subsets and their progeny retain the capacity to generate all other subsets and re-establish the parental EML cell 
population, while maintaining their distinct phenotypic and functional properties, (a) FACS purified parental CD34 ne8 , CD34'°™, CD34 lnt and CD34 hl8h 
EML cell subsets (cell purity > 99% and viability > 95%) were cultured for 14 days (n = 3). (b) Flow cytometry analysis of CD34 expression on day 14 
progeny of each parental EML cell subset. Day 14 progeny of each parental cell subset was sorted into separate sets of "primary" CD34 neg , CD34 low , CD34"" 
and CD34 h ' 8h subsets (n = 3). (c) The overall differentiation pattern of each of the primary cell subsets was analyzed using CFC asays. Shown here as an 
example is the differentiation pattern of the CD34 ne8 , CD34 low , CD34' nt and CD34 h ' 8h progeny of the parental CD34 high cell subset (n = 15). Data are shown 
as mean ± SD. (d) Sets of primary CD34 ne8 , CD34'°™, CD34 mt and CD34 hl8h cell subsets, generated by each parental cell subset, were again cultured 
independently for another 14 days, (e) Representative flow cytometry analysis of CD34 expression profiles on day 14 progeny of each primary EML CD34 
cell subset has shown that all primary cell subsets re-established the CD34 expression profile of the parental EML cell population. Shown here as an 
example is CD34 expression on the progeny of primary CD34 ne8 , CD34 low , CD34 int and CD34 l, ' 8h cells generated by parental CD34 low cell subset (n = 3). 



expression on day 14 total progeny of each parental EML CD34 cell 
subset has again shown that each subset generated all other cell 
subsets and re-established the parental EML cell population 
(Fig. 7b). Day 14 progeny of each parental EML cell subset was sorted 
into separate sets of "primary" CD34 nee , CD34 low , CD34 int and 
CD34 hlBh subsets (99% purity) (Fig. 7b). CFC assays have shown that 
the overall differentiation patterns of each of these primary cell sub- 
sets resembled the differentiation patterns of parental cell subsets 
(Fig. 7c). Each of purified primary EML CD34 cell subsets was cul- 
tured independently for 14 days (Fig. 7d). Analysis of CD34 express- 
ion profile on their progeny showed that each of the primary EML 
CD34 cell subsets again generated all other CD34 cell subsets and 



re-established the broad expression range of CD34 marker typical for 
the parental EML cell population (Fig. 7e). 

Cumulatively, these studies revealed that on a population level all 
EML CD34 cell subsets exhibit cell-autonomous interconvertible 
properties, where each subset self-renews and generates all other 
subsets ultimately re-establishing the complete parental population 
of EML cells. Furthermore, the interconversion capacity is main- 
tained by all serially generated EML cell subsets through multiple 
cell generations. 

Analysis of the interconversion capacity of CD34 neg , CD34 low , 
CD34 in, and CD34 hl8h single cell clones. Although our data clearly 
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show that the CD34 neg , CD34 low , CD34 ,nt and CD34 h ' Bh EML cell 
subsets contain cells with interconversion capacity, the question 
remained whether all or only a fraction of cells in each CD34 
subset have that capacity. 

The instructions for maintenance of EML cell line (obtained 
directly from Dr. Tsai) 30 and our own long-term experience show 
that EML cells grow optimally when maintained at a concentration of 
—2 X 10 5 cells/ml. Culturing EML cells below the concentration of 
0.5 X 10 5 cells/ml slows down the cell growth, increases apoptosis 
and leads to population drift and culture collapse. When cultured as 
single cells in 96-well plates only a fraction of EML cells (<20%) 
generated some progeny (—10 to <500 cells) which rarely expanded 
to the level necessary for proper flow cytometry analysis. Similarly, 
only a small percentage of total EML cells and purified CD34 neB , 
CD34 1 "", CD34 int and CD34 high EML cells forms single cell-derived 
undifferentiated blast colonies (CFU-Blast) 30 when plated at a low 
cell density (10 3 cells/ml) into methylcellulose cultures supplemen- 
ted with recombinant SCF (Supplementary Fig. S8). Although these 
colonies can be extracted and cultured in EML media, the majority of 
extracted CFU-Blast colonies do not expand sufficiently in liquid 
cultures. 

In an attempt to analyze interconvertibility of EML CD34 cell 
subsets at a single cell level, we undertook a different experimental 
approach. The CD34 neg , CD34 low , CD34 inl and CD34 hi s h cell subsets 
were purified from GFP-expressing EML cells and single GFP + cells 
from each subset were cultured together with wild type GFP~ EML 
cells as "support" cells in a mixed culture system. 

GFP + EML cells were generated by transduction of wt EML cells 
with murine stem cell virus (MSCV) -based bicistronic retroviral 
vector, with IRES driving expression of GFP. EML cells were trans- 
duced at several different MOI (transduction efficiency of 70-85%), 
and 72 hours after transduction, GFP + EML cells were purified by 
sorting (Supplementary Fig. S9a). Multiple GFP + EML cell clones 
were generated, expanded and early passages cryopreserved. All 
GFP + EML cell clones have shown long-term (>6 months), stable 
expression of GFP, and displayed the same growth characteristics as 
wt EML cells (data not shown). 

Flow cytometry analysis did not reveal any significant differences 
in the expression pattern of c-kit, Sca-1, Flk-2 and CD34 markers 
between wt EML cells and GFP + EML cells, (Supplementary Fig. 
S9b,c and Fig. la). The analysis of the cell cycle status with PI has 
reproducibly shown that the percentage of wt EML cells and GFP + 
EML cells in G0/G1 and S/G2/M phases was almost identical 
(Supplementary Fig. S9d and Fig. 2b). In the hematopoietic CFC 
assays wt EML and GFP + EML cells displayed almost identical multi- 
lineage differentiation capacity (Supplementary Fig. S9e). These 
results indicate that transduction with MSCV and GFP expression 
had no effects on the HSC marker expression pattern, cell cycle 
profile and differentiation capacity of GFP + EML cells. Moreover, 
GFP + EML cells show interconvertible properties identical to wt 
EML cells (data not shown). 

The CD34 neg , CD34 low , CD34 ,nl and CD34 hish cell subsets were 
purified from GFP+ EML cells by cell sorting (Fig. 8a). The high 
purity cells from each subset were deposited at a concentration of 
£1 cell/well into individual round bottom wells within 96-well 
plates; with each well containing 1 X 10 4 wt GFP~ EML cells in 
100 ul of EML media (Fig. 8a). The cells in each well were expanded 
by passaging the cell progeny sequentially into 24-well plates, 6-well 
plates, T25 culture flasks and finally into T75 culture flasks. 

Cell aliquots from each of these mixed cultures were analyzed by 
flow cytometry to identify all cultures containing GFP + progeny 
from single CD34 ne s, CD34 low , CD34 ,nt and CD34 hieh cells (Fig. 8a). 
In these experiments -21% of single GFP + CD34 nee cells, -28% of 
single GFP + CD34 low cells, -48% of single GFP + CD34 tot cells and 
—57% of single GFP + CD34 hlgh cells generated various levels of pro- 
geny (Fig. 8b). 



GFP+ progeny of each CD34 neg , CD34 low , CD34 int and CD34 h ' gh 
single cell clone was purified from mixed cultures by cell sorting and 
cultured independently for 14 days to expand the cells for flow cyto- 
metry analysis (Fig. 8a). While the progeny of the CD34 mt and 
CD34 h ' eh single cell clones expanded very well, the progeny of the 
CD34"' e and CD34 low single cell clones did not expand efficiently to 
generate enough cells for proper flow cytometry analysis. 

Nevertheless, flow analysis of the progeny of multiple CD34 lnt and 
CD34 h ' eh single cell clones has shown that each clone generated all 
other CD34 cell subsets and re-established the broad expression 
pattern of CD34 marker typical for complete parental EML cell 
population (Fig. 8c,d and data not shown). 

Discussion 

To gain further insight into the heterogeneity and interconvertibility 
of multipotent hematopoietic cells, we analyzed in more depth the 
heterogeneity and interconvertible properties of cells within the mul- 
tipotent cell line EML. The expression pattern of a panel of HSC 
markers revealed an extensive phenotypic heterogeneity among 
EML cells, which resembles heterogeneous expression of various 
markers on HSCs and their subsets 6,15,42-49 . Using broad range of 
expression of the CD34 marker EML cells can be reproducibly 
divided into CD34 ne e, CD34 low , CD34 int and CD34 hlgh cell subsets. 
Remarkably, all EML CD34 cell subsets are phenotypically quite 
heterogeneous as well, and display distinct expression pattern of 
various HSC markers (Sca-1, CD27, CD150, CD201, VEGFR1). In 
that regard, each EML CD34 cell subset can be further divided into 
distinct subpopulations based on the co-expression patterns of vari- 
ous HSC markers. 

Functional characterization of EML CD34 cell subsets has revealed 
that they display distinct multilineage differentiation capacity, cell 
cycle profiles, CFSE labeling profile and proliferation kinetics when 
compared against each other. Moreover, the immunophenotypic and 
functional blueprints of each EML cell subset suggest a closer asso- 
ciation between the CD34" e6 and CD34 low cell subsets on one end and 
the CD34 tot and CD34 hiBh cell subsets on the other end. 

By establishing phenotypic and functional signatures of EML 
CD34 cell subsets our studies have revealed even more extensive 
cellular heterogeneity and complexity within EML cell line, and sup- 
port the notion that EML cell line consists of multiple distinct cell 
subsets. 

The microarray analysis and comparison of hematopoietic tran- 
scriptional signatures in CD34" eB EML cells and combined CD34 ml 
and CD34 hlgh subsets did not reveal significant differences in the 
expression of the majority of (1) genes important for overall hema- 
topoiesis and HSC maintenance and function, (2) hematopoietic cell 
surface markers, receptors and adhesion molecules, and (3) genes 
associated with erythroid, megakaryocyte, myeloid and lymphoid 
lineage specification and development 50 52 . These results are very 
similar to findings of earlier studies, which reported that expression 
levels of most genes are quite similar in Sca-1~CD34~ and Sca- 
1 + CD34 + EML cells 34,41 . 

Previously, Sca-l low EML cells were found to express significantly 
higher levels of GATA-1 and lower levels of PU.l transcript and 
protein than Sca-l h ' Bh EML cells, although at the single cell level 
the Sca-l low EML cells are quite heterogeneous for GATA-1 express- 
ion 32,34,41 . Division of Sca-l low EML cells into CD34~ and CD34 + 
subsets revealed that <80% of Scal low CD34" EML cells and less 
than 5% of Scal low CD34+ cells express GATA-1 34 . 

Our studies have shown that expression of GATA- 1 transcript is 
the highest in the CD34 nee and CD34'°™ cell subsets and significantly 
lower in the CD34" 11 and CD34 h ' eh subsets. In contrast, the expression 
of PU. 1 transcript is the highest in CD34 h ' eh subset and is significantly 
lower in the CD34 1 "', CD34 low and CD34 neg subsets. 

Western analysis has revealed that the expression of GATA-1 
protein is higher in the CD34 neg and CD34 low cells, while the express- 
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Figure 8 | Analysis of interconversion capacity of CD34 ne8 , CD34 low , CD34"" and CD34 hlgh single cell clones, (a) Experimental approach to analyze 
interconversion capacity of GFP + CD34 ne6 , CD34 low , CD34" 11 and CD34 hleh single cell clones using mixed culture system with wild type GFP~ EML cells as 
"support" cells, (b) Frequency of CD34 1 " 6 , CD34 low , CD34 1 " 1 and CD34 hl8h single cells that generated GFP + cell progeny in mixed culture system, (c) 
Representative flow cytometry density plots showing the CD34 expression pattern on parental GFP+ EML cells and the progeny of two independently 
generated CD34 lnt and CD34 h ' gh single cell clones, showing re-establishment of the broad expression range of CD34 marker typical for parental EML cell 
population (n = 3). (d) A representative overlay of the expression range of CD34 marker on GFP+ parental EML cells and the progeny of two different 
subclones of independently generated CD34'" 1 and CD34 hlgh single cell clones (n = 6). 



ion of PU.l protein is the highest in the CD34 hlBh and CD34 tot cells. 
However, the differences in the expression of GATA-1 protein 
between EML CD34 cell subsets were not as prominent as the ones 
observed between Sca-l low , Sca-l med and Sca-l M8h EML cells, and Sca- 
,iow C D34~ and Seal 1 ™ CD34 + EML cells 32 ' 34 . 

That dissimilarity is most likely due to differences in the cellular 
composition of Sca-l low , Sca-l med , Sca-l MBh , Sca-l low CD34~ and 
Scal low CD34 + EML cells on one end and EML CD34 cell subsets 
on the other. For example, the CD34 neg cell subset encompasses the 
cells displaying the full range of Sca-1 expression (i.e. Sca-l 1 " 8 , Sca- 
l'°", Sca-l med and Sca-l hi8h cells) (Supplementary Fig. S4), and con- 
tains previously analyzed Sca-l low , Sea- 1 1 ™ 1 , Sca-l hlgh , Sca-l low CD34" 
and Scal low CD34 + EML cells. Thus, the level of GATA-1 expression 
found in CD34 neB cell population is an average of different levels of 



GATA-1 expression found previously in various Sca-1 EML cell 
subsets. 

Transcription factor C/EBPa controls developmental fate of mul- 
tipotent hematopoietic cells by promoting myeloid and inhibiting 
erythroid differentiation 55,56 . Both earlier and our own microarray 
studies indicated that the expression of C/EBPa transcript is (1) 
higher in Scal + CD34 + EML cells than in Scar CD34~ EML cells 41 , 
and (2) higher in CD34 h,gh and CD34" 1 ' cells than in CD34 neg and 
CD34 low cells (Fig. 3). Although our qRT-PCR analysis revealed 
highly significant difference in the expression of C/EBPa transcript 
between CD34 high and CD34 in, cells on one end and CD34 ne8 and 
CD34 low cells on the other end, there was no significant difference 
in the expression of C/EBPa protein between all EML CD34 cell 
subsets (Fig. 3). 
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The observed differences between the expression levels of GATA- 
1 and C/EBPtx transcripts and proteins in EML CD34 subsets could 
be the consequence of an unequal rate of GATA-1 and C/EBPcx 
protein translation or differential rate of protein degradation in these 
cell subsets. 

Cumulatively, these studies have discovered differential express- 
ion of some of the key transcriptional regulators of lineage specifica- 
tion and commitment between EML CD34 cell subsets, and the 
grouping of CD34 ne 8 and CD34 low cells versus CD34 hi < ih and 0034"" 
cells in terms of their expresion pattern of lineage-associated tran- 
scription factors. 

Moreover, these results reaffirm the notion that EML CD34 cell 
subsets are heterogeneous and encompass cells in different stages of 
lineage specification and commitment. In that regard the CD34 nt ' s 
EML cell subset is functionally heterogeneous and consists of several 
cell subpopulations, including a subpopulation of cells that express 
high levels of GATA-1, are more erythroid committed but exhibit 
low differentiation rate into erythroid and myeloid lineages. 

On the other end, the expression pattern of lineage-associated 
transcription factors and more efficient diffrentiation into myeloid 
and erythroid lineages suggest that the CD34 mt and CD34 h ' eh subsets 
contain cells that are either oligopotent and give rise to multiple 
lineages in response to cytokines or are committed only to erythroid 
or myeloid lineages. 

Using phenotypic and functional signatures of EML CD34 cell 
subsets as a platform, we interrogated their maintenance during 
multiple cell generations. Temporal analysis of the primary progeny 
of each purified EML CD34 cell subset established that as a cell 
population each CD34 subset exhibits cell-autonomous intercon- 
vertible properties, and contains cells that both self-renew and gen- 
erate other cell subsets, ultimately re-establishing complete parental 
EML cell population (Supplementary Fig. SlOa). 

Further testing of interconvertible characteristics of each EML 
CD34 cell subset through serial cell sorting, culture and analysis of 
their successive progeny over several generations demonstrated that 
all sequentially generated CD34 cell subsets maintain their intercon- 
vertible properties and the capacity to self-renew, generate all other 
CD34 subsets and re-establish the parental EML cell population, 
while maintaining cell subset-specific phenotypic and functional 
features (Supplementary Fig. SlOb). These findings support the view 
that EML cell line consists of multiple interconvertible cell subsets 
which co-exist in a dynamic but stable equilibrium through multiple 
cell generations. 

Equally important, perpetual maintenance of interconvertible 
EML cell subsets in the absence of cellular and molecular compo- 
nents of BM HSC niches 59 (with an exception of SCF), infers that 
similar to intrinsically regul ated interconvertibility of ES cell sub- 
sets 1 3 ' 20,21 , the dynamic equilibrium of interconvertible EML cell 
subsets is homeostatically maintained and regulated at least in part 
by currently unknown cell-intrinsic mechanisms. 

While all EML CD34 cell subsets clearly contain cells with inter- 
conversion capacity, the question remained whether all or only a 
fraction of cells in each subset have that capacity. 

Chang et al. used methylcellulose cloning method to generate 
random single cell-derived clones from total EML cells that 
expressed varying levels of Sca-1. Over a period of 8 weeks these 
clones regenerated the original broad expression profile of Sca-1, 
typical for EML cell line 32 . Pina et at sorted single cells from Sca- 
1'°, Seal 10 CD34~ and Seal 1 " CD34+ subsets into 96-well plates and 
analyzed their so called culture-reconstituting potential over a period 
of several days. In that study 17% of Sca-1 10 , 2.5% of Seal 10 CD34" 
and 35% of Scal'° CD34 + cells generated more than 100 cells, which 
supposedly also reconstituted broad Sca-1 expression 34 . 

Using a mixed culture system we generated multiple CD34 neB , 
CD34 low , CD34 int and CD34 Mgh single cell clones. However, we were 
able to sufficiently expand only the progeny of CD34 mt and CD34 h ' 6h 



single cell clones, which exhibited interconversion capacity by gen- 
erating all other EML CD34 cell subsets and re-establishing the broad 
expression pattern of CD34 marker typical for complete parental 
EML cell population. 

Our studies are clearly showing that on a cell population level both 
CD34"' e and CD34 low cell subsets exhibit interconvertible properties. 
However, the culture of purified CD34 nt ' 6 and CD34 low cells has con- 
sistently shown their delayed kinetics of regenerating the parental 
EML cell population. Moreover, that kinetics in mixed culture with 
other cell subsets being present remains to be determined. Thus, 
although the progeny of CD34 nt ' 8 and CD34 low single cell clones did 
not expand sufficiently to determine their capacity to regenerate 
parental EML cell population, that outcome does not mean that 
single cells from these subsets do not possess interconvertible 
properties. 

Collectively, our findings further highlight the phenotypic and 
functional heterogeneity of EML cell subsets, and bring additional 
insight into the characteristics and extent of EML cell interconvert- 
ibility. In addition, together with previous studies, our results further 
emphasize the biological relevance of EML cell line as a unique 
experimental model to study the process of interconversion and cell- 
intrinsic and extrinsic pathways that regulate interconversion among 
multipotent hematopoietic cells. 

While the previous and our studies demonstrated interconversion 
capacity of EML cells at a single cell level to a varying extent, the 
actual identity and phenotypic and functional features of these cells 
remain to be fully elucidated. 

Purification of specific cell subpopulations with defined immuno- 
phenotypes within each EML CD34 cell subset now allows identifica- 
tion of exact cell fractions in each subset that have interconversion 
capacity. In that regard, the cellular barcoding 60 of phenotypically 
defined EML cell populations with documented interconversion 
capacity could enable more precise in vitro tracking of the progeny 
of single cells, and help to more accurately analyze the process of 
interconversion in a mixed culture system with multiple subsets 
present at the same time. This approach could also provide a plat- 
form for elucidating intrinsic molecular mechanisms regulating the 
process of interconversion between different states. 

Better defined heterogeneity of EML cell subsets, purification of 
specific cell subpopulations based on their immunophenotypes, and 
the molecular analysis on a single cell level 15 could provide a high 
resolution insight into the developmental and lineage commitment 
stages, and developmental hierarchy of cells composing each EML 
cell subset. 

Methods 

Cell lines and cell culture. The Stem Cell Factor (SCF)-dependent EML cell line was 
kindly provided by Dr. Steve Collins (Fred Hutchinson Cancer Research Center, 
Seattle, WA). EML cells were maintained in IMDM with L-glutamine (Gibco), and 
supplemented with 20% (vol/vol) equine serum (Hyclone), 1% Penicillin/ 
Streptomycin (Sigma) and 10% (vol/vol) conditioned media from the BHK/MKL cell 
line (constitutively expressing recombinant soluble SCF) as previously described 31,37 . 
EML cells were maintained at a density of ~2 X 10 5 cells/ml, and were passaged every 
other day. The BHK/MKL cell line was maintained in DMEM with L-glutamine 
(Gibco), 10% fetal bovine serum (Hyclone), and 1% Penicillin/Streptomycin (Sigma). 
Conditioned media was collected from confluent BHK/MKL cell cultures at 48 and 
96 hours, centrifuged at 1200 rpm for 10 minutes, filtered through a 0.45 um 
membrane (Nalgene) and stored at — 20°C. 

Flow cytometry analysis and cell sorting. EML cells were harvested, washed with 
phosphate buffered saline (PBS), and incubated with antibodies at concentrations of 
1-2 ug per 10 5 cells in 100 ul of PBS for 40 min on ice, followed by double wash in 
PBS. Centrifugation steps were carried out at 1200 rpm for 10 min at 4 L 'C. The cells 
were analyzed live or after fixation in 10% neutral buffered formalin. The following 
antibodies were used: IgG2a and IgG2b isotype controls (eBioscience), anti-c-kit APC 
(eBioscience, clone 2B8), anti-Sca-1 PE-Cy7 (eBioscience, clone D7), anti-Flk-2 PE 
(eBioscience, clone A2F10), anti-CD49b PE-Cy7 (eBioscience, clone DX5), anti- 
CD34 FITC (eBioscience, clone RAM34), anti-CD34 eFluor660 (eBioscience, clone 
RAM34), anti-CD150 APC (eBioscience, clone 9D1), anti-CD201 AlexaFluor647 
(eBioscience, clone mRCR-16), anti-VEGFRl APC (R&D Systems, clone 141522) 
and anti-CD27 PE-Cy7 (eBioscience, clone LG.7F9). Before labeling of cells with 



SCIENTIFIC REPORTS | 4:5199 | DOI: 1 0. 1 038/srep05 1 99 



13 



multiple antibodies, different fluorochromes were checked for compatibility and 
compensation utilizing the BD Fluorescence Spectrum Viewer. The expression of 
various markers on cells was analyzed on BD LSR II or BD LSRFortessa flow 
cytometers (BD Biosciences) using the live cell gate. Flow cytometry data were 
analyzed with Flow Jo software (Tree Star). EML cells and cell subsets were sorted on a 
FACSAria II sorter (BD Biosciences). Only the cells from sorts yielding the highest 
cell purity (>99%) and viability (>95%) were used for experiments. 

Cell cycle analysis and proliferation kinetics. The cell cycle profile of total EML cells 
and purified CD34 neg , CD34 low , CD34 int and CD34 high EML cell subsets was analyzed 
by Propidium Iodide (PI), Hoechst 33342, and Ki-67/7-AAD staining. For staining 
with PI (BIOTIUM) the cells were resuspended in cold 70% ethanol, washed in PBS, 
treated with 0.1 mg/ml RNase A, and incubated with 0.1 p.g/ul PI. For staining with 
Hoechst 33342 (SIGMA) the cells were resuspended in cold PBS/1% formaldehyde, 
and incubated with 5 ug/ml Hoechst at 37°C for 15 minutes. For analysis of Ki-67 
expression, the EML cells were first labeled with anti-CD34 FITC antibody. Staining 
of labeled EML cells with Ki-67-eFluor450 {eBioscience, clone B56) followed the 
eBioscience Protocol B for Staining Intracellular Antigens for Flow Cytometry (http:// 
www.ebioscience.com/media/pdf/best-protocols/staining-intracellular-antigens- 
for-flow-cytometry.pdf), utilizing the Foxp 3 /Transcription factor staining buffer set 
(eBioscience), followed by addition of 1 : 100 ul of 7-AAD (BioLegend) in PBS. 

To examine their proliferation kinetics, total EML cells and EML cell subsets were 
labeled with CFSE (Carboxyfluorescein diacetate succinimidyl ester) (Invitrogen) by 
adding 1.5 uM CFSE in 1 ml of 0.1% FBS/1X PBS for 8 minutes at room temperature, 
followed by addition of 1 ml of warm FBS and 10 minute incubation in 37°C water 
bath. The cells were washed 3 times in 2% FBS/1X PBS and placed in culture. The 
division of CFSE labeled EML cells was evaluated at 24 hour time points for 5 days. 
Aliquots of CFSE labeled cells at each time point were incubated with anti-CD34 
eFluor660 antibody (eBioscience) as described. At each time point only the viable 
cells, as determined by staining with PI or LIVE/DEAD® Fixable Yellow Dead Cell 
Stain Kit (Invitrogen), were analyzed. 

Hematopoietic colony-forming unit assays. To assess their multilineage 
differentiation capacity, total EML cells and purified CD34 neg , CD34 low , 0034"' and 
CD34 Mgh EML cell subsets were plated in 2 ml of 0.3% low melting temperature 
agarose (SeaPlaque, Cambrex) with 2X IMDM (Gibco), 20% heat inactivated equine 
serum (Hyclone) and 10% SCF- conditioned medium in replicate 6-well plates 
(10 3 cells/ well) 37 ' 38 . The media was supplemented with: 0.1 ng/ml of recombinant 
mouse interleukin 3 (IL-3) (R&D Systems) for generation of granulocyte- 
macrophage colony forming units (CFU-GM); 8 U/ml of recombinant human 
erythropoietin (Epo) (Ortho Biotech) for generation of burst-forming erythroid units 
(BFU-E); and 3 ng/ml of recombinant mouse Thrombopoietin (Tpo) (R&D Systems) 
for generation of megakaryocyte colony forming units (CFU-Meg). Cultures were 
maintained at 37°C, 5% C0 2 for 7-10 days. The BFU-E, CFU-GM and CFU-Meg 
colonies were counted 7-10 days after plating. Data from multiple independent 
experiments (n — 15) are shown as mean ± standard deviation (SD). 

Microarray expression analysis. To examine molecular differences between EML 
cell subsets we analyzed the transcriptomes in three separately purified samples of 
CD34 neg cells and combined CD34 int and CD34 high subsets, refrerred to in the text as 
CD34 + cells. The microarray expression analysis was performed by Ocean Ridge 
Biosciences (ORB, Palm Beach Gardens, FL) using Agilent mouse microarrays. The 
mRNA samples were prepared from total EML cells and purified CD34 neg and CD34 + 
EML cells using Micro -FastTrack™ 2.0 mRNA and mRNA DIRECT Isolation Kits 
(Life Technologies). Ocean Ridge Biosciences performed additional purification of 
mRNA, sample labeling, microarray hybridization and data analysis. The results were 
provided by ORB in Log2 transformed form and are presented as mean ± s.d. (n — 3). 
The <1. 5-fold difference in expression was not considered significant. To identify 
differences in the hematopoietic transcriptomes between EML CD34 cell subsets we 
specifically examined the expression of (1) genes important for overall hematopoiesis 
and HSC maintenance and function, (2) hematopoietic cell surface markers, 
receptors and adhesion molecules, and (3) genes associated with erythroid, 
megakaryocytic, myeloid and lymphoid lineage specification and development. 

Quantitative real-time reverse transcription PCR analysis. mRNA was isolated 
using the Micro -FastTrack™ 2.0 mRNA and mRNA DIRECT Isolation Kit (Life 
Technologies). Equal amount of mRNA from each cell sample was reverse 
transcribed using High Capacity cDNA Reverse Transcription Kit (Applied 
Biosystems). Quantitative RT-PCR reactions (qRT-PCR) reactions were performed 
in an Agilent Technologies Stratagene Mx3005P for 40 cycles (95 °C for 30 seconds, 
57-62' C for 1 minute, and 72°C for 1 minute). The specificity of PCR products was 
verified by melting curve analysis and by running qRT-PCR products on 2% agarose 
(SeaKem, Cambrex) gels and analyzing them on Alphalmager 3400 (Alpha 
Innotech). The results were normalized for expression of HPRT as a control, and are 
shown as fold difference compared with the levels (set to 1) detected in the sample 
with the lowest expression. The sequences of primers used in qRT-PCR analysis are as 
follows: GATA-1 (sense- CAGGGCAGAATCCACAAACT, antisense- 
TCCTCTGCATCAACAAGCC); PU.l (sense- TGACTACTACTCCTTCGTGG, 
antisense- GATAAGGGAAGCACATCCGG); C/EBPoc (sense- 
C A AG AG CCG AG AT AAAGC, antisense- GTCATTGTCACTGGTCAA); JAM4 
(sense-CTGTTGCTGTTGTGCCTC, antisense- TTGCCTC ATCCGAACTGT A) ; 



GAPDH (sense- TGCACCACCAACTGCTTAGC, antisense- GGCATGGACT- 
GTGGTCATGAG). 

Western blotting. The cells were lysed in 100 ul of modified RIPA buffer (150 inM 
NaCl, 1% NP-40, 150 mM Tris, pH 7.5, 1 mM EDTA), supplemented with protease 
inhibitor (Roche), PhosStop (Roche), 10 mM HEPES and 200 mM PMSF. 25-40 jig 
of protein from each sample was run on 12% SDS-PAGE gels (Bio-Rad Mini Protean- 
TGX) and transferred to nitrocellulose membranes (Bio-Rad Trans-Blot). Prior to 
transfer the gels were stained with Coomassie Brilliant Blue (Bio-Rad). Rat anti- 
GATA-1 (clone N6), rabbit anti-PU.l (clone T-21), goat anti-JAM4 (clone M-13), 
goat anti-rat IgG-HRP and goat anti-rabbit IgG-HRP antibodies were purchased 
from Santa Cruz Biotechnology. Rabbit anti-C/EBPa monoclonal antibody (clone 
EP709Y) was purchased from Epitomics. Membranes were incubated overnight with 
1 : 200 anti-GATA-1, 1 : 500 anti-PU.l, 1 : 500 anti-C/EBPa and 1 : 500 anti-JAM4 
antibody in 5% BSA at 4°C, washed and then incubated overnight at 4°C with 
1 : 10,000 goat anti-rat IgG-HRP or goat anti-rabbit IgG-HRP secondary antibodies. 
Proteins were detected with TMB Membrane Peroxidase Substrate System (KPL) for 
detection of HRP. Human monocytic cell line U-937 (ATCC, CRL-1593.2) were used 
as positive controls for expression of C/EBPa protein. Membranes were stripped and 
incubated with rabbit anti-GAPDH antibody (Sigma- Aldrich). The levels of GATA- 
1, PUT, C/EBPa, JAM4 and GAPDH proteins in each sample were measured by 
densitometric scanning and analyzed with ImageJ software. 

Statistical analysis. Statistical significance was determined using a two-tailed 
unpaired Student t-test (Prism Software, Graphpad). Data are expressed as mean ± 
standard deviation (SD). The number of independent samples and experiments is 
provided in figure legends. Statistical significance and P values are indicated by 
asterisks and specified in figure legends. P values of <0.05 were considered 
statistically significant. 
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